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High dynamic range imaging by pupil single-mode filtering and 
remapping 
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ABSTRACT 

Because of atmospheric turbulence, obtaining high angular resolution images with a high 
dynamic range is difficult even in the near infrared domain of wavelengths. We propose a 
novel technique to overcome this issue. The fundamental idea is to apply techniques devel- 
oped for long baseline interferometry to the case of a single-aperture telescope. The pupil 
of the telescope is broken down into coherent sub-apertures each feeding a single-mode 
fiber. A remapping of the exit pupil allows interfering all sub-apertures non-redundantly. A 
diffraction-limited image with very high dynamic range is reconstructed from the fringe pat- 
tern analysis with aperture synthesis techniques, free of speckle noise. The performances of 
the technique are demonstrated with simulations in the visible range with an 8 meter telescope. 
Raw dynamic ranges of 1 : 10 6 can be obtained in only a few tens of seconds of integration time 
for bright objects. 

Key words: atmospheric effects - turbulence - instrumentation: adaptive optics - techniques: 
high angular resolution - techniques: interferometric - (stars:) planetary systems 



1 INTRODUCTION 

Most astronomical sources have a small angular extension. The 
study of their spatial intensity distribution thus requires a high an- 
gular resolution. Among these, exoplanets require in addition very 
high dynamic ranges to be directly imaged near their host star. High 
angular resolution and high dynamic range are the key. Because of 
diffraction, the angular resolution of a light collector linearly in- 
creases with its diameter. Larger telescopes therefore provide larger 
angular resolutions. However, phase aberrations due to atmospheric 
turbulence prevent from reaching the diffraction limit. The adap- 
tive optics (AO) technique, in which a deformable mirror corrects 
the corrugated wav efront, allows to rest ore the diffraction limit in 
spite of turbulence iRousset et al Jl 1990T> . Although an elegant and 
effective solution, it provides a limited sky coverage in absence of 
a polychromatic laser guide star and is still difficult to apply to 
the next generation of very large telescopes at short wavelengths. 
Post-processing techniques such as speckle imaging and aperture 
masking, allow restoring diffraction limited images from short ex- 
posure i mages in which turbulent fluctuations are frozen. Speckle 
imaging <Labevridll97Ch has a limited dynamic range and aper- 
ture masking iHaniff etalll987l) only makes use of a small fraction 
of the pupil. Here we propose to combine the advantages of these 
two techniques and to spatially filter the beams to achieve high dy- 
namic range diffraction limited images even with extremely large 
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telescopes at optical and near-infrared wavelengths. We first dis- 
cuss the issue of imaging through the turbulent atmosphere in the 
optical in Section|2| We present our instrument concept in Section|3| 
and discuss the anticipated performances in Sectionfl] Conclusions 
are drawn in Section|5| 



2 IMAGING AT OPTICAL WAVELENGTHS THROUGH 
ATMOSPHERIC TURBULENCE 

The theory of imaging establishes that a telescope acts as a low- 
pass filter on an image, the optical transfer function (OTF) being a 
decreasing function from the zero (spatial) frequency down to the 
telescope cut-off frequency D/X with D the pupil diameter and A 
the wavelength (Figure 0. In presence of turbulence and for long 
exposures, the cut-off frequency is reduced to ro/A when ro < D 
hence a spatial resolution limited to A/rrj. The Fried parameter, ro, 
is the diameter of a coherent patch over which the variance of phase 
is 1 rad 2 . AO partly restores the OTF up to D/ X as long as the in- 
coming wavefront can be analyzed on a time scale shorter than the 
coherence time of turbulent phase. In practice, this requirement is 
fulfilled at near-infrared wavelengths as the coherence time and ro 
increase with wavelength. Correction of turbulent phase is there- 
fore more demanding in the visible and no large astronomical tele- 
scope has been equipped with a visible AO system yet, despite the 
potentially higher angular resolution. Correction of turbulent phase 
for future extremely large telescopes at near-infrared wavelengths 
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Figure 1. 8 m telescope Optical Transfer Functions. Grey continuous line: perfect telescope without turbulence. All three other OTFs are for excellent 
turbulence conditions (seeing of 0.4 arc second). Black continuous line: OTF with the system proposed in this paper. The width of the bins is proportional to 
the diameter of a sub-pupil. The OTF in each bin is constant and is equal to the weighted average of the filtered optical transfer function. Black dotted line: 
with adaptive optics correction. Black dashed line: with the speckle interferometry technique. Left: at a visible wavelength. Right: in the near infrared. The 
quality of the OTF improves with wavelength making performance at short wavelengths quite poor. With the pupil remapping system, the OTF is mechanically 
constant making the telescope in turbulent conditions a quasi-perfect telescope thus restoring diffraction limit and high dynamic range. 



will also be more difficult as the number of coherent patches will 
be more than an order of magnitudes larger. With the speckle tech- 
nique, the OTF is also restored up to the same limit D/X. In both 
cases however, the higher the spatial frequency, the larger the en- 
ergy loss. This is due to the larger decorrelation of turbulent phases 
between the more distant parts of the pupil. 

OTFs obtained in different conditions with an 8 m telescope 
have been simulated and are presented in Figure Q Simulations 
were made in the visible (A = 0.55 /^m) and in the near infrared 
(A = 1.65 /im). For turbulent conditions, an excellent seeing of 
0.4 arcsecond has been chosen. The OTF of a perfect telescope (no 
turbulence and no aberrations) is compared to the OTFs obtained 
with adaptive optics, with the speckle technique and with the sys- 
tem proposed in this paper (the width of the bins is the width of the 
sub-pupils). The performance of adaptive optics and speckle in- 
terferometry are highly turbulence-dependent and therefore highly 
time and wavelength-dependent. In practice, turbulent conditions 
can change over short time scales making the calibration of the 
OTF difficult. The situation improves with wavelength making per- 
formance at short wavelengths quite poor. With the pupil remap- 
ping system, the OTF is mechanically constant making the tele- 
scope in turbulent conditions a quasi-perfect telescope thus restor- 
ing diffraction limit and high dynamic range (see Section|3j. 

The losses of coherence are responsible for speckle noise in 
the image that limits the dynamic range with both speckle imaging 
(to about 1 : 100) and AO making it difficult to reach 1 :10 5 dynamic 
ranges required for exoplanetary system imaging <Chelliil2004 
ICavarroc et alj|200r3) . They are also responsible of a smaller con- 
trast in the image on smaller spatial scales. In the classical imaging 
approach, the pupil is massively redundant as a same vector can be 
drawn between several pairs of points in the pupil (a given spatial 
frequency is measured different times in the pupil). In the OTF, a 
frequency component is the sum of vectors of complex amplitude 
produced by different pairs whose phases are the phase differences 
between two points of the incoming wavefront. The degradation of 
the OTF is therefore induced by these incoherent additions in a re- 
dundant pupil. The technique of non-redundant aperture masking 



overcomes this issue: a mask is placed in a pupil plane to select a 
set of sub-pupils organized as a non-redundant array. In the instan- 
taneous OTF, each discrete spatial frequency is therefore generated 
by a single sub-aperture pair and there is no coherence loss when 
the sub-pupils are point-like or far smaller than r . As a conse- 
quence, the OTF is a set of peaks of unit amplitude beside the zero 
frequency peak. An image can then be reconstructed through the 
closure phase technique. This technique has been demonstrated on 
the sky on several telescopes and spectacular astronomical results 
have been obt ained such as the imaging and di scovery of the pin- 
wheel nebula <Tuthill. Monnier & Danchill999l) . 



3 PUPIL REMAPPING AND SPATIAL FILTERING WITH 
SINGLE-MODE FIBERS 

The shortcoming of non-redundant aperture masking is that only a 
small amount of the pupil can be used thus leading to a very poor 
sensitivity. The diameter of the sub-pupils therefore needs to be in- 
creased as much as possible to about tq. This however reduces the 
average coherent energy per frequency bin to about 36% and the 
OTF is fluctuating and difficult to calibrate. This yields a degraded 
dynamic range in the reconstructed image. This is a well-known is- 
sue in long baseline interferometry. It has been solved with single- 
mode fibers. Their fundamental property is that the spatial distri- 
bution of the lightwave phase is by nature flat across a waveguide 
section and the transported lightwave is the refore b y definition per- 
fectly coherent. The FLUOR instrument ICoude du Foresto et alj 
[l998j), based on this technology, has demonstrated that the coher- 
ent energy in long-baseline interferometry is fully restored leading 
to more accurate visibility measurements and to potentially higher 
dynamic range reconstructed images. The key is that the phase cor- 
rugations over each pupil of the interferometer are traded against 
intensity fluctuations which can be monitored and calibrated at the 
pace of turbulence. The same principle can be applied to aperture 
masking. The lightwave from each sub-pupil feeds a single-mode 
fiber. The net effect is that the frequency peaks in the OTF of a 
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Figure 2. Instrument set-up. The input pupil is mapped with a lenslet array. Each lenslet focuses a sub-pupil beam into a single-mode fiber. In the output pupil, 
the fibers are placed at the focus of another lenslet array. Fibers are distributed in a non-redundant array. Each pair of sub-pupils in the input pupil yields a 
single spatial frequency component in the output pupil. Spatial frequency components with different turbulent phases are therefore not mixed thus preserving 
optimum imaging performances. 



non redundant array are automatically of unit amplitude whatever 
the size of the sub-apertures. The u se of fibers in non-redun dant 
aperture masking has been proposed IChang & BuscheJl998h . 

We propose to further take advantage of the use of fibers to 
rearrange the distribution of sub-apertures in the exit pupil of the 
instrument (Figure[2j. The input pupil is mapped with a lenslet ar- 
ray. Each lenslet focuses a sub-pupil beam into a single-mode fiber. 
The focal length of lenslets is chosen for optimum flux injection, 
the maximum coupling efficiency is 78% for a step-index fiber. The 
full entrance pupil of the telescope can now be used and turned into 
a non-redundant array in the exit pupil of the instrument. In the out- 
put pupil, the fibers are placed at the focus of another lenslet array 
and recast into a non-redundant array. Each pair of sub-pupils in the 
entrance pupil of the telescope is associated to a single sub-pupil 
vector in the exit pupil. The different phases of redundant pairs are 
therefore not mixed by the interferometric combination and the full 
coherent energy is preserved. The OTF modulus in the exit pupil 
is equal to 1 for each exit spatial frequency or, equivalently, the 
OTF in the entrance pupil is maximum for each bin (Figure^. As 
far as the source is concerned, the measurement of the modulus of 
each spatial frequency component is unbiased by speckle noise and 
limited by ph oton and detector no ise only. The rotation shearing in- 
terferometer iMariottiet all 19921) and the coding of non-redundant 
sub-a rrays with differential wavefront tips and tilts iMonnier et alj 
120041) are an alternative to achieve non-redundant pupil recombina- 
tions. However, both setups do not include spatial filtering and the 
achievable dynamic range is therefore limited. 



4 PERFORMANCES OF THE TECHNIQUE 

Although modal filtering by fibers imposes a flat wavefront in each 
sub-pupil, the average differential phases between sub-apertures 
are preserved and the exit wavefront is a step function. The clo- 
sure phase technique allows eliminating these phase errors as the 
sum of phases over a triangle of sub-apertures is not affected by 
any sub-aperture phase error. It provides a fraction of 1 — 2 /TV (N 
is the number of sub-apertures) of the phase information. In our 
case however, the high redundancy of the entrance pupil is useful. 



Closure phases measured in the exit pupil yield the above fraction 
of phase information. Under high entrance redundancy, the num- 
ber of exit closure phases is larger than the number of indepen- 
dent phases therefore providing the full phase information. This 
can be understood with a system as simple as N = 6 hexagonal- 
shaped sub-pupils arranged on a hexagon. The number of pairs is 
•Npairs = N(N — l)/2 = 15 and the number of independent 
triangles is iVtriangics = (N - 1)(N - 2)/2 = 10 . However, 
most pairs form a same vector and the number of different vectors 
is equal to 9. In fact, the number of measured closure phases is 
larger than the amount of available object phase information. This 
effect increases with the number of sub-apertures in such a com- 
pact array. Thanks to th e redundancy, the flux coupled in each fiber 
can also be measured lLacour. Thiebaut & Perrin|l2006l) . Eventu- 
ally, our system allows measuring all complex spatial frequency 
components of the object as sampled by the array of sub-pupils. 
The accuracy on moduli and phases can be directly derived from 
the achievement of single-mode long baseline interferometers: a 
fraction of a percent for moduli and better than a degree for phases. 
Better performances can be expected with a single-telescope exper- 
iment. Using conservative SV = 1% and 5<f> = 1°, an estimate of 
^the dy namic range is obtained from the formula iBaldwin & Hanifu 

DR=J - pa "' s (1) 

V {SV/Vf +5<f> 2 

For an unresolved object (V = 1), the dynamic range is DR m 
50 X Npairs- Assuming a 132 sub-pupil array, a dynamic range 
of 1:4600 is achieved per snapshot. Since speckle noise has been 
eliminated, different snapshots are uncorrelated and the dynamic 
range for S snapshots is equal to 1:4.6 x 10 3 x y/S . With 10 4 
snapshots, a dynamic range of 1:4.6 x 10 l> is obtained in the recon- 
structed image. The relevance of this formula is illustrated by Fig- 
ure 3. The observation of a multi-component object under good see- 
ing conditions (0.6 arc seconds) with a 132-element pupil remap- 
ping system at the focus of an 8 meter telescope at 630 nm wave- 
length (bandwidth = 60 nm) has been simulated. No adaptive op- 
tics correction is applied and only 6% of the flux is coupled in 
the fibers in average. Performances in injected flux would at least 
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Figure 3. Simulated observation and image reconstruction. Simulations of the observation of a multi-component object under good seeing conditions (0.6 arc 
seconds) with a 132-element pupil remapping system at the focus of an 8 meter telescope at 630 nm wavelength (bandwidth = 60 nm). No adaptive optics 
correction is applied. The object is a star surrounded by a disk whose integrated flux is a hundredth of that of the star. A decreasing exponential law has been 
chosen for the disk flux. Two planets orbit the star and are dimmer by factors of 10 3 (Planet A) and i0 4 (Planet B) with respect to the central object. The 
closest planet is located 6 diffraction limits away from the star. The star image is saturated and the diffraction limit can be assessed on an unresolved planet 
image. The magnitude of the star is respectively V = 0, 5, 10 from right to left. Each image is the sum of 10000 snapshots of 4 ms amounting to a total exposure 
time of 40 s. Images are recorded in the photon noise limited regime with a photon counting camera. The dynamic ranges directly measured in the images are 
of 1:10 4 , 1 :10 s and 1:10 6 as predicted. The faintest component is clearly detected in the V=5 and images. The speckle noise has been eliminated by the 
instrument and detection is photon noise limited. 



double with an upstream near-infrared optimized adaptive optics 
correction. The object is a star surrounded by a disk whose inte- 
grated flux is a hundredth of that of the star. A decreasing expo- 
nential law has been chosen for the disk flux that is only a 1000 th 
of the star brightness at a distance equal to three times the diffrac- 
tion limit. Two planets orbit the star and are dimmer by factors of 
10 3 (Planet A) and 10 4 (Planet B) with respect to the central ob- 
ject. The closest planet (Planet A) is located 6 diffraction limits 
away from the star. The star image is saturated in the displayed 
images and the diffraction limit can be assessed on an unresolved 
planet image. The magnitude of the star decreases from the left 
to the right by increments of 5 in magnitude (a factor of 100 in 
brightness, stars from the left to the right have visual magnitudes 
of 10, 5 and 0, respectively). Each image is the sum of 10000 snap- 
shots of 4 ms amounting to a total exposure time of 40 s. The 
frame rate has been chosen so that the piston is frozen during an 
acquisition to prevent the fringe pattern to be blurred and to keep 
the fringe contrast larger than 99% (fringe contrast accuracy better 
than 1%). Using the parameters of the simulations, the maximum 
allowed frame time is 12 ms. Should the frame rate be smaller than 
100 Hz then the final dynamic range could be smaller than those 
achieved here. Techniques to measure the contrast loss due to pis- 
ton could possibly be developed to relax this constraint and allow 
larger individual exposure times. Images are recorded in the pho- 
ton noise limited regime with a photon counting camera. Images 
ar e reconstructed from the interfero metric data with the algorithm 
of lLacour. Thiebaut & Perrir] l200(j|) . The dynamic ranges directly 
measured in the reconstructed images are of 1:10 4 , 1:10 s and 1:10 6 
as predicted (thanks to twice more accurate visibilities than in the 
above calculation). Th e faintest component is clearly de tected in 
the V=5 and images. lLacour. Thiebaut & Perrinl l2006t) demon- 
strate that speckle noise is eliminated and detection is photon and 
detector noise limited as anticipated when using single-mode fibers 
and measuring the differential phases between sub-pupils. This is 
where the gain in dynamic range comes from. The dynamic range 



can therefore still be increased by accumulating more images to 
reduce the impact of the noise. 

These performances do not take into account systematics that 
could impair the actual performance of the instrument. As long as 
the bandwidth is small enough (i.e. the coherence length of the ra- 
diation is larger than the rms piston) the instrument is immune to 
bandwidth (chromatic) effects, piston and higher order turbulence 
effects. Problems may arise from imperfections in the realization 
of the optics and in the alignment of the instrument such as polar- 
ization and dispersion whose effects will reduce the fringe contrast 
and alter the fringe phase. However such static effects can be mit- 
igated by calibrating the data with data acquired on a calibrator 
star unresolved by the instrument. First experiments will allow to 
study and discover those systematics and measure their impact on 
the technique. 

Since the image is reconstructed from the measurement of spa- 
tial frequency components, the noise is white in the final image. 
This is a drawback with respect to direct imaging techniques for 
which the source photon noise is essentially concentrated in the 
central star image and can be rejected thanks to a coronograph. A 
nulling technique can however be used by applying a phase pattern 
to the fiber network to cancel light wherever needed. The corono- 
graph is however of no use for multiple bright sources and for ex- 
tended objects whereas our technique can be applied to any astro- 
nomical source. The advantage of our technique is the perfection 
of the OTF hence the absence of speckle noise. A better detectivity, 
or equivalently a higher dynamic range, could be achieved in dif- 
ferent ways. The point spread function can be accurately modeled 
and subtracted from the reconstructed images and advanced faint 
object algorithms could be used to detect planets. Alternatively, the 
visibility data can be directly fitted with a model of the object with- 
out reconstructing an image to search for faint planets around stars 
to avoid reconstruction noise (the spatial frequency plane coverage 
is not perfect). 
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The high dynamic range can be as well obtained on more ex- 
tended objects. The instantaneous field of view is primarily limited 
by two effects: the lobe of each single-mode fiber matched to X/d 
where d is the diameter of a sub-pupil ; and the interferometric 
field of view. The latter comes from the non-homothetic relation 
between the entrance and exit pupils of the instrument. In usual 
non-Fizeau combination schemes (i.e. for which the exit and en- 
trance pupils are not homothetic), the exit pupil is densified with 
respect to the entrance pupil to reduce the number of peaks (or 
fringes) in the fringe pattern to concentrate the coherent flux in the 
central peak of width X/D with D the size of the entrance pupil. 
The most extrem e case of such a set-up is the hyper-telescope of 
lLabevrieiiT99r3) . In this case, the interferometric field of view is a 
fraction 7 of X/d with 7 the densification factor. In the present in- 
strument, the pupil is not densified but diluted by a large amount to 
reach a non-redundant configuration. The exit diffraction pattern is 
composed of peaks whose characteristic size is A/_B max with B max 
the longest baseline in the exit pupil. Those peaks are spread over 
a single sub-aperture diffraction pattern of width X/d 1 . The inter- 
ferometric field of view then depends upon spectral bandwidth. As 
a matter of fact, the wave carrier of the fringe pattern has a width 
at least equal to that of the diffraction pattern of a sub-pupil. For an 
off-axis object (off-axis distance a) to contribute to the fringe pat- 
tern of the on-axis source, it is necessary that the coherence length 
of the radiation in baseline units (i.e. divided by the maximum base- 
line D) be as large as a. Hence the classical formula for the field 
of view of a co-axial beamcombiner: 

a ^mx- (2) 

The field of view is maximum (i.e. the interferometric field of view 
matches the fiber lobe) when: 

AA=(^)xA. (3) 

For the simulation presented in this paper, this holds if the spectral 
bandwidth is roughly equal to 60 nm hence our choice. The field 
of view can therefore be adjusted according to the characteristics 
of the astrophysical targets. If needed, the field of view can still be 
increased by a 2-D mapping of the sky if a field of view larger than 
the fiber lobe is necessary. 

The use of an upstream near-infrared optimized AO is very 
advantageous as it allows to increase the size of each sub-pupil 
and therefore the sensitivity. The single-mode aperture masking 
imager can therefore be used on the largest telescopes at visible 
wavelengths behind a near-infrared optimized AO or on future ex- 
tremely large telescopes either at visible or infrared wavelengths 
and will be able to provide very high angular resolution and dy- 
namic range images. Last but not least, since the wavefront after 
filtering with fibers is a step function, the imager could also ac- 
tively control wavefront errors, the actuators being fiber stretchers. 
In this active mode, the instrument could be a simpler AO version 
for future extremely large telescopes. 



5 CONCLUSION 

We have proposed a novel technique to obtain high angular res- 
olution and high dynamic range images through a turbulent at- 
mosphere whatever the wavelength, be it in the visible or in the 
infrared domains. The telescope pupil is broken into sub-pupils 
each injected into single-mode fibers to perfectly spatially filter the 
beams corrugated by atmospheric turbulence. A non-redundant out- 
put pupil is formed to avoid mixing phases from different parts of 
the input pupil. We have shown that the optical transfer function 
of a telescope could be restored quasi perfectly with this technique 
and that speckle noise is eliminated. Raw dynamic ranges of 1:10 s 
can be obtained on bright sources in a few tens of seconds of inte- 
gration time. This makes the technique very promising for the di- 
rect imaging of exoplanets around nearby stars but also for the use 
of current ground based telescopes and future extremely large tele- 
scopes at very high angular resolution even at short wavelengths. 
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1 This can be easily demonstrated e.g. in the case of a compact square array 
of square sub-apertures. When applying a magnification factor to the grid 
spacing of the array of sub-aperture centers, peaks appear for factors larger 
than 1 whereas there is a single peak for a magnification factor of 1 (i.e. 
no magnification). This can be generalized to a non-regular array of sub- 
apertures as is the case for a 2-D non-redundant array but is beyond the 
scope of this paper. An example of such diffraction pattern can be found in 
Figure 3 o flTuthill et. all |20QO|) . 



